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Abstract Models for Majorana neutrino masses can be classified according to the level in 
perturbation theory at which the effective dimension five operator LLHH is realized. The 
possibilities range from the tree-level up to the three-loop level realizations. We discuss 
some general aspects of this approach and speculate about a model independent classifica- 
tion of the possible cases. Among all the realizations, those in which the effective operator 
is induced by radiative corrections open the possibility for lepton number violation near -or 
at- the electroweak scale. We discuss some phenomenological aspects of two generic ra- 
diative realizations: the Babu-Zee model and supersymmetric models with bilinear R-parity 
violation. 
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1 Introduction 

Neutrino experiments have firmly demonstrated neutrinos are massive and have non- vanishing 
mixing angles among the different generations 11]. Since in the standard model neutrinos are 
massless this experimental results are a clear evidence of beyond standard model physics. 
From a general point of view Majorana neutrino masses can be generated by adding to the 
standard model Lagrangian the non-renormalizable effective Lagrangian [2J 

ifg ^ ^ ^5 = _L ^iT^H) {HiT2L) , (1) 
Anp Anp 

where L and H are the lepton and Higgs electroweak doublets. The presence of the ef- 
fective Lagrangian in ([T]l guarantees that after electroweak symmetry breaking neutrinos 
acquire Majorana masses. Extensions of the standard model in which neutrinos have Majo- 
rana masses realize the effective operator 0'^ in different ways. However, among all these 
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possibilities there are subsets that have a common feature, namely the order in perturbation 
theory at which the operator is realized. As will be discussed in sec.[2]a rather general clas- 
sification of neutrino mass models can be done by gathering together these different subsets. 

Making an exhaustive list of possible realizations of ff^ and their corresponding phe- 
nomenology is far out of the scope of the present discussion. Thus, instead of using this 
approach we will discuss two classes of generic models: models in which neutrino masses 
are generated radiatively at the 2-loop level and supersymmetric models with broken R- 
parity. In the former case we will stick to the Babu-Zee model t3ji4J, while for the latter we 
will discuss bilinear R-parity breaking models with a neutralino LSP (Si- 



2 The ff's operator and its different realizations 

At Anp the heavy physical degrees of freedom, that when integrated out yield the dimension 
five {d = 5) effective operator ff^, are no longer decoupled. Once specified they define a 
particular model for Majorana neutrino mass generation, so different models lead to differ- 
ent realizations of 6^ . As pointed out in the introduction these realizations can be classified 
according to the order in perturbation theory at which ^5 is generated, namely tree-level 
(^1""), one-loop (^5""'), two-loops (^5^^) and three-loops {ff'i^^'), where (.. denotes the 
number of loops in each case. Three-loop level realizations require order one Yukawa cou- 
plings implying that models based on ^g-"* are, in general, not consistent with neutrino data 
once the requirement of perturbativity of the corresponding couplings is imposed. 

A general and model independent classification of the tree-level ^g"" and one-loop 
level ^1^' realizations of ^5 has been carried out in ref. f6l and such an analysis could 
be, in principle, extended to the two and three loop order realizations. Regardless of the 
order at which the operator arises the procedure is based on the determination of all the 
possible gauge invariant renormalizable vertices within the loop for all the different possible 
topologies. In this procedure the only gauge quantum numbers that are fixed are those of the 
external legs {LLHH) while those of the physical degrees of freedom flowing in the loop 
are free and should be fixed by the requirement of gauge invariance. Another approach for 
the same sort of classification is by using all the possible AL = 2 effective operators up to 
certain d. This method has been implemented in ref. \1 1 for effective operators up to d = 11. 
The difference between these two approaches is that while the former focus exclusively on 
the 1^5 operator the later covers this operator but also higher dimensional operators as e.g. 
^7 ~ (LLHH){H'^ H). Note that these higher dimensional effective operators {d > 5) can 
give a dominant contribution to neutrino masses only if the leading effective operator ff^ is 
forbiden due to a new symmetry. These cases have been throughout analysed in ref. 1 8 1 and 
we have nothing more to add here. 

The tree-level realizations of ^5 correspond to the different type of seesaw models (type- 
I f9), type-II 1 10] and type-Ill 1 1 1 ]) whereas for the 1-loop cases the number of possibilities 
(models) is much more larger (is determined by the different SU (3) x SU (2) x [/( 1 ) assign- 
ments of the internal degrees of freedom flowing in the loop). Examples of include 
the Zee model \\T\, models with scalar leptoquarks 113] and models with extra scalars and 
fermions with nontrivial color charges |14J- Examples of two-loop realizations include ex- 
tended scalar sectors as in the case of the the Babu-Zee model |3 ,4] and models with scalar 
leptoquarks fTlJlSl. Up to our knowledge three-loop level realizations rely on extensions of 
both the scalar and fermion sectors, examples can be found in references f7] |16l . 

In models in which the operator ^5 is generated beyond the tree-level the lepton number 
breaking scale can readily be at or around the electroweak scale. Thus, models embedded in 
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Fig. 1 Feynman diagram for the 2-loop Majorana neutrino masses in the model of 13 141 . 



such realizations usually lead to testable predictions in either high-energy or high-intensity 
experiments. In what follows we will discuss two main cases, the Babu-Zee model and 
supersymmetric bilinear R-parity violating models, paying special attention to some of their 
phenomenological implications. 



3 Two loop realization: the Babu-Zee model 

In this model the standard model scalar sector is extended by the addition of two new scalars, 
/i+ and both singlets under SU{2). Their couplings to standard model leptons is given 
by 

^ = /a/J {l^^aLCl^'^d^ah^ + ^a^ (eaRCep«)A:++ + h.c. (2) 

Here, L/, are the standard model (left-handed) lepton doublets, eg the charged lepton singlets, 
a,jS are generation indices and e,j is the completely antisymmetric tensor. Note that / is 
antisymmetric, while is symmetric. Assigning L = 2 to hr and eq. ^ conserves 
lepton number. Lepton number violation in the model resides only in the following term in 
the scalar potential 

^ = -\lh"^h^k— (3) 

Here, ^ is a parameter with dimension of mass. 

The setup of eq. (|2]l and eq. ((S} generates Majorana neutrino masses via the two-loop 
diagram shown in fig. ([TJ. The resulting neutrino mass matrix can be expressed as 



with summation over x,y implied. The parameters a)„ are defined as a)cy = mjixyniy, with 
OTv the mass of the charged lepton (v Following 1 17] we have rewritten h^a = h'^a 
hap = 2h'^p. J'ir) finally is a dimensionless two-loop integral given by 



- fdxf 

Jo Jo 



1 y(l-v) 
dy—, log ■ (5) 



x+{r—l)y + y^ x + ry 



For non-zero values of r, J^{r) can be solved only numerically. We note that for the range of 
interest, say 10^^ <r < ■^{r) varies quite smoothly between (roughly) 3 < ^{r) < 0.2. 
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Fig. 2 Conservative lower limit on the branching ratio Br(fi — > ey) as a function of the charged scalar mass 
mi, for normal hierarchy (left plot) and inverted hierarachy (right plot). The three lines are for the current 
solar angle sin^ G12 best fit value (full fine) and 3 fj lower (dashed line) and upper (dot-dashed line) bounds. 
Other parameters fixed at sin^ 023 = 0.5, sin^ 613 = 0.040 and Am\^^^ = 2.0 ■ 10"'' eV^. 

3.1 Flavour violating charged lepton decays 

Phenomenological tests of this model have been studied in I17II18I[T91 . Among all of them 
those involving jU — > ey can be regarded as the most stringent ones. In ref. IIS! it has been 
shown that the corresponding decay branching ratio for this process can be written as 

with e = fell flit and m/, the mass of the singly charged scalar. Figure [2] shows the resulting 
lower limit on Br{^ — )• ey) as a function of m/, for the case of normal and inverted hierar- 
chies. Note that the horizontal solid line indicates the upper limit set by the MEGA experi- 
ment [201 and not the new one placed by the MEG experiment, Br{ji — )• ey) < 2.4 x 10"'^ 
at 90% C.L. [21 1- Using the updated limits the constraints on the singly charged scalar mass 
would be even more stringent that the ones quoted here. 

In summary, in this model Br{^ — > ey) > 10"'-^ is guaranteed for singly charged scalar 
masses smaller than 590 GeV (5.04 TeV) for normal (inverse) hierarchical neutrino masses, 
and larger or even much larger branching ratios are expected in general. Thus, a non- 
observation of this process in the next few years, at least for the case of inverse hierarchy, 
would certainly remove most of the motivation to study this model. 

4 Bilinear R-parity violating supersymmetry 

Bilinear R-parity violation (BRpV) provides an intrinsically supersymmetric framework for 
Majorana neutrino masses (for a review see 1221 '). In these models the superpotential in- 
cludes, in addition to the MSSM terms, also the term 

H'BRpv = e,i,4. (7) 

This term breaks not only R-parity but also lepton number in all three generations. In order 
to have a consistent model a soft SUSY breaking term has to be added to the scalar potential. 
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namely 

^BRpV = BiEiLiHu . (8) 

The presence of these terms induce a non-vanishing vacuum expectation value for the sneu- 
trinos (v; = (v,)) that give rise to a mixing among neutral gauginos and Higgsinos with 
neutrinos. Due to this mixing one of the neutrinos acquires mass. The effective neutrino 
mass matrix reads 

4det(M^o) ^ ' 

where A,- = jixv,- + v^e (with /i the Higgsino mass term and = {Hd}) and M^o is the neu- 
tralino mass matrix. The other two neutrinos acquire mass from one-loop corrections involv- 
ing W, Z and scalar loops, being the bottom-sbottom and tau-stau loops the most important 
contributions f23l. Thus, the BRpV model is an example of a model in which neutrino 
masses arise from different realizations of the ^5 operator. 

The atmospheric and reactor angles are approximately given by the rotation angles that 
diagonalize the tree level mass matrix my"' : 

tan2e23«^, tan^ei3« ,2^! ,2 ■ (^O) 

The solar angle instead is obtained once the one-loop corrections are taken into account. A 
remarkable feature of the bilinear R-parity breaking model is that the same parameters that 
determine neutrino physics also control the decay patterns of the LSP, thus it is always pos- 
sible to establish a set of correlations between the LSP decay branching ratios and neutrino 
observables. These correlations can be used as a tool to know whether BRpV is responsible 
for the origin of neutrino masses. 



4.1 LSP decays and neutrino observables 

Once R-parity is broken the LSP is unstable and decays to standard model fermions. Thus 
astrophysical contraints on its nature do not hold any more and in principle any supersym- 
metric particle can be the LSP. The decay patterns of all possible LSPs and their relations 
with neutrino observables in the context of bilinear R-parity breaking models have been 
analysed in f5l!241. Additional studies in more generic models including trilinear R-parity 
violating couplings have been carried out for slepton and sneutrino LSPs |25,26]. Here we 
will highlight some of the main features of this "program" in bilinear R-parity breaking 
models assuming the lightest neutralino to be the LSP. The following discussion is based on 
ref. la. 

The presence of the bilinear R-parity violating parameters induce not only a mixing 
between neutralinos and neutrinos but also a mixing between charginos and charged leptons, 
charged Higgses and sleptons, CP-even (CP-odd) components of the neutral Higgses and 
the corresponding CP-even (CP-odd) components of the sneutrinos (23"]. All together these 
mixings determine the three-body leptonic, semi-leptonic and invisible neutralino decays: 
Xi Vilpk' X't ^N^' Xi Viqq and Xi ^ V/V^Vj. 

Given the above discussion it should now become clear that once the bilinear R-parity 
breaking parameters are constrained by neutrino experimental data the lightetest neutralino 
decays are expected to be constrained as well. In fact, it turns out that the constraints imposed 
by neutrino physics assure: (i) the neutralino always decay inside the detector; (ii) the decay 
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Fig. 3 Ratios of decay brandling ratios for semi-leptonic neutralino decays and their correlations with the 
atmospheric and reactor angles (S)- 



branching ratio for neutralino invisible decays never exceeds 10%; (iii) different ratios of 
decay branching ratios are strongly correlated with neutrino mixing angles. Figure [3] shows 
the corresponding correlations for the semi-leptonic final states fiqq', Tcjq' and eqq' . From 
these results and the measured neutrino mixing angles it can be established that if BRpV is 
the mechanism responsible for the origin of neutrino masses and the lightest neutralino turn 
out the be the LSP the following measurements should be expected at LHC: 

BR{x'-^^fiqq')-mXi^^qq') and BR(x['^e^9')«BR(Z?^Mw')- (H) 

In summary, in bilinear R-parity breaking models the decay patterns of the LSP are 
strongly correlated with neutrino mixing angles. These correlations allow to set constraints 
on the different decay branching ratios of the LSP that in turn can be used to prove whether 
these models are responsible for the origin of neutrino masses and mixings. 



5 Conclusions 

From a general perspective Majorana neutrino masses can be accounted for by the effective 
dimension five operator ^g. We have argued that the different "incarnations" of this oper- 
ator can be classified according to the order in perturbation theory at which the operator is 
realized. In principle by using group theoretical arguments one could make, in a model inde- 
pendent way, an exhaustive list of all the possibilities at each order and up to the three-loop 
levelQ. 

Models in which ffs is realized radiatively (ff^^^^, where £ denotes the numer of loops) 
rely on TeV scale physics. Thus, an obvious question is whether this new physics, and there- 
fore the origin of neutrino masses, can be proved at e.g. the LHC. As illustrative examples 
we have discussed what we consider two benchmark models: The Babu-Zee model and su- 
persymmetry with bilinear broken R-parity. 
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